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Abstract 

We study charmless B + meson decays to the pA.TT + n~ final state using a 605 fb -1 data sample 
collected at the Y(4«S) resonance with the Belle detector at the KEKB asymmetric-energy e + e~ 
collider. There are significant signals found with the pA mass peaking near threshold. The observed 
branching fraction for non-resonant B + — > pA7r + 7r~ is (5.92±§;!f (stat.) ± 0.69(syst.)) x 1(T 6 with 
a significance of 9.1 standard deviations. We also observe the intermediate three-body decay 
B + -> pAp° with a branching fraction of (4.78±° %l{stat.) ± 0.60(syst.)) x 1CT 6 and a significance 
of 9.5 standard deviations, and find a hint of a B + — > pA/2(1270) signal. No other intermediate 
three-body decay is found in this study. 

PACS numbers: 13.25.Hw, 13.30 Eg, 14.40.Nd 
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I. INTRODUCTION 



The large data samples accumulated at the B factories allows us to study the properties 
of baryonic B meson decays, which are less well understood than B decays into mesons. An 



barvomc 



interesting hierarchy has been established experimentally in branching fractions of 
B decays, namely B(B° -> pA c ~7r+7r) > B(B + -> pA c ~7r+) > B(B° -> pA~) Q, 3, S, El|. 
The above decays proceed via b — > c tree diagrams. The corresponding hierarchy has not 
yet been confirmed in charmless baryonic -B decays, which presumably proceed via b — > s 
penguin or b — > it tree diagrams. So far, many three-body charmless baryonic B decays 
have been observed but only very stringent upper limits have been set for quasi-two-body 
decays One intriguing experimental finding in charmless baryonic three-body decays is 



that all the baryon-antibaryon mass distri 



Dutions peak near threshold. Many theoretical 



investigations, using QCD counting rules or pole models in the factorization 

approach, have been carried out for charmless baryonic B decays with pA in the final state. 
These theoretical results on the branching fractions and the shape of the baryon-antibaryon 
mass spectrum agree reasonably well with the experimental findings. However, the proton 



angular distribution of the threshold peak observed in B° —> pAir~ 10[ seems to violate 
the short- distance b — > s picture and awaits a theoretical explanation. Since the A hyperon 
could be a useful tool to probe the b — > s process jf], [nj and the related three-body decays 
can be used for a T- violation study [ljj, it is of general interest to search for additional 
charmless baryonic decay modes. 

In this paper, we present for the first time a study of the four-body charmless baryonic 
decay B + — > pAn + ir~ [k]] in order to test the multi-body hierarchy in the charmless case. To 
investigate the threshold enhancement effect found in three-body decays, we study partial 
branching fractions as a function of the baryon-antibaryon mass. We also search for possible 
intermediate three-body decays to the same pAir + ir~ final states, such as pAp°. 

II. EVENT SELECTION AND RECONSTRUCTION 
A. Data Samples and the Belle Detector 

The data sample used in this study corresponds to an integrated luminosity of 605 fb _1 
and contains 657 xlO 6 BB pairs collected with the Belle detector on the T(4S) resonance 
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at the KEKB asymmetric-energy e + e~ (3.5 GeV and 8 GeV) collider [jjj]. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter composed of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux-return 
located outside of the coil is instrumented to detect K° L mesons and to identify muons. The 



Belle detector is described in detail elsewhere 



15|- 



B. Selection Criteria and B Meson Reconstruction 

The selection criteria for the final state particles of B + — > pAtt + 7t~ are based on in- 
formation obtained from the tracking system (SVD and CDC) and the hadron identifica- 
tion system (CDC, ACC, and TOF). We use the same requirements as those described in 



Ref. [10[. Candidate B mesons are identified with the following two kinematic variables in 
the center-of-mass (CM) frame: the beam-energy-constrained mass M bc = \/E^ eam — p B , 
and the energy difference AE = Eb — -Ebeam, where -Ebeam is the beam energy, and ps and 
Eb are the momentum and energy, respectively, of the reconstructed B meson. We define 
a candidate region as 5.20 GeV/c 2 < M bc < 5.3 GeV/c 2 and -0.1 GeV < AE < 0.5 GeV. 
The lower bound in AE is chosen to exclude all backgrounds due to baryonic B decays 



with higher multiplicities. From a GEANT 



Monte Carlo (MC) simulation, the signal 



peaks in a signal box defined by the requirements 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and 
\AE\ < 0.05 GeV. 

C. Background Suppression 

After the above selection requirements, the background in the candidate region arises 
predominantly from continuum e + e~ — > qq (q — u, d, s, c) processes. We suppress jet-like 
continuum background relative to the more spherical BB signal using a Fisher discrim- 
inant 17] that combines seventeen event shape variables as described in Ref. We 
then optimize the coefficients separately in seven different missing-mass regions defined in 
Ref. [l9j to improve the signal-to-background ratio. Probability density functions (PDFs) 
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for the Fisher discriminant and the cosine of the angle between the B flight direction and the 
beam direction in the T(45 l ) rest frame are combined to form the signal (background) likeli- 
hood L s (£b). The signal PDFs are determined using signal MC simulation; the background 
PDFs are obtained from the sideband region of the data: 5.20 GeV/c 2 < M\ x < 5.26 GeV/c 2 
or \AE\ > 0.1 GeV. We require the likelihood ratio 1Z = C S /(C S + £&) to be greater than 
0.85 for all seven missing-mass regions. This value is determined by optimizing n s /y/n s + rib 
as a function of TZ, where n s and n& denote the expected numbers of signal and background 
events in the signal box, respectively. We assume a signal branching fraction of 10~ 5 to 
estimate n s and use the data sideband events to determine n&. 

To ensure that the decay process is genuinely charmless, we apply a charm veto. 
B + — > pAii + 7i~ candidate events with 2.10 GeV/c 2 < M A7r + < 2.32 GeV/c 2 are removed to 
avoid background from B + — > pA c ir + and B + — > pS c °, S c ° — > A c ir + with A c — > Air~ 
decay. From MC simulation, there are events from B° — > pAir~ in the candidate region. For 
simplicity, we remove B candidates if the corresponding reconstructed Mb c and AE, using 
only the three daughters (i.e., pAir~), are in the signal box. This selection removed about 
4.8% of candidate events. The contribution of the B background component with S — > A7 
has a AE distribution that is different from signal.. This is included in the systematic un- 
certainty from PDF modeling by comparing the fit results with and without this background 
component in the fit. If there are multiple B candidates in a single event, we select the one 
with the best x 2 value for the pir + 7r~ vertex fit. The fraction of multiple B candidate events 
is 18.8%. The systematic errors due to multiple B candidates are described later. 



III. EXTRACTION OF SIGNAL 

In order to obtain the signal yield for the pAn + ix~ final state, we perform an unbinned 
extended likelihood fit that maximizes the likelihood function 

L= M 

N 

X n(^pA7r+7r--PpA7r+7r- + ^iqPqq), 
i=l 

where is the number of total events. A r p ^ 7r+7r _ and N qq - are fit parameters representing 
the numbers of pAir + 7T~ signal events and continuum background events, respectively. Each 
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function P is a PDF expressed as uncorrelated products of shapes of the M^ c and AE 
distributions: P = P]\i hc x Pae- 

For the PDFs related to B decays, we use a Gaussian function to represent -Pju bc an d a 
double Gaussian for Pae with parameters determined from a MC signal simulation of B + — > 
pAn + ir~ phase space. To model the PM bc continuum background, we use a parameterization 
that was first employed by the ARGUS collaboration, f{M\ )C ) oc x\/l — x 2 e~^ l ~ x2 \ where 
x is M bc /_E bcam and £ is a fit parameter j^. The Pae continuum background shape is 
modeled by a normalized second-order polynomial whose coefficients are fit parameters. 




FIG. 1: B signal yields obtained as a function of the pA mass. The yields are corrected for the 
mass-dependent efficiency. The solid curve shows a fit with a threshold function. 



We study mass spectra for the pA, pix~ , p7r + , Att~ , Att + , 7r + 7r~, pAn~ , pA7r + , 
piT + 7r~ and A7r + 7T~ combinations using the B signal yields obtained as functions of those 
masses. There is a clear enhancement near threshold in the pA mass for signal candi- 
dates. We fit the pA mass distribution with a threshold function: / t h r oc (AM j) s x 

e [ Cl x(AM pX ) + c 2 x(AM pX f+c 3 x(AM p7 ,f]^ where ^ M __ = M ^_ _ ^ _ ^ ^ ^ ^ ^ ^ fit 

parameters. The result is shown in Fig. [IJ The excess in the region 3.5 GeV/c 2 < M p ^ < 4.0 
GeV/c 2 will be investigated in a study of systematic effects. The only other resonance- like 
structure found is in the n + ir~ mass spectrum where we observe a clear p° signal. We cor- 
rect the fitted B signal yields by the efficiencies to obtain the total B(B + — > pA7r + 7r~)Tot = 



7 



(11. 28lo?2 ^ 1-03) x 10~ 6 , where non-resonant and possible intermediate resonance decays 
are all included. 

In order to investigate three-body decays such as B + — > pAp° and B + —>■ pA/ 2 (1270), we 
perform an unbinned extended likelihood fit that maximizes the likelihood function 

e -( N P A. + n~ + N P A P + N P Af 2 + N w) 

L = M 

N 

X n(^pA7r+7r-^pA7r+7r- + ^pApPpAp 
i=l 

+^pA/2^pA/2 + NqqPq-q), 

where N is number of total events. Here N p -^ n+w -, N p -^ p0 , N p -^^ 1270 ^, and N g g are fit 
parameters representing the yields of non-resonant pA7r + 7r~, pAp° and pA/ 2 (1270) sig- 
nal and continuum background contributions , respectively. Each function P is a PDF 
expressed as uncorrelated products of shapes of the M bc , AE and M n + n - distributions: 

P = Pmu c x Pae x Pm + _ • 

For the PDFs related to B decays, we use a Gaussian function to represent PM hc and 
a double Gaussian for Pa_b with parameters determined from MC signal simulation. For 
the continuum background PDF, we use an ARGUS function to model the Pu hc distribu- 
tion while the Pae distribution is modeled by a normalized second-order polynomial whose 
coefficients are fit parameters. 

The PDFs for the M n + n - mass distributions in B + — > pAp° and B + — > pA/ 2 (1270) are 
smoothed histograms obtained from MC simulation using a sample with the pA shape fixed 
from data and Breit-Wigner forms for the p° or f 2 resonances. The subsequent decay angular 
distribution is assumed to be flat. The PDF Pm , _ for B + — > pAtt + tt~ is also a smoothed 
histogram from a three-body phase space MC sample, which again follows a threshold shape 
in the pA mass spectrum for data. The continuum background PDF for M^^- is modeled by 
the sum of smoothed histogram functions taken from Breit-Wigner distributions for p° and 
fo resonances, and a threshold function: Pm + - = r\X P p + r\X Pf + (1 — n — r 2 ) x P t hr and 
P thr oc (AM 7r+7r -) s x e [ ClX ( AM -+--) +C2X ( AM -+--) 2 ] where AM n+n - = M 7r+7r --2m 7r ; n, r 2 , s, a 
and c 2 are fit parameters. 
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TABLE I: Systematic errors (in %). 



Source 


pA7T + 7T 


pAp° 


pA/ 2 (1270) 


Tracking 


6.0 


5.7 


5.7 


Proton ID 


4.2 


4.2 


4.2 


Charged Pion ID 


1.3 


1.7 


1.4 


A reconstruction 


3.3 


2.8 


3.0 


B{A -» P TT+) 


0.8 


0.8 


0.8 


B{ P /h - 7T+7T') 


- 


0.2 


2.8 


1Z Selection 


2.4 


2.4 


2.4 


Multiple Counting 


3.0 


3.0 


3.0 


PDF Modeling 


7.3 


9.0 


9.6 


Number of BB Pairs 


1.4 


1.4 


1.4 


Total 


11.7 


12.6 


13.4 



IV. RESULTS 

Figure [2] shows the fit results for B + — > pAtt + tt~. The resulting signal yields are 
I67.8+2I?, 131.2±^|;| and 39.11^ with statistical significances of 9.1, 9.5, and 3.0 for 
B + — > pA7r + Tc~, B + — > pAp° and B + — > pA/ 2 (1270), respectively. The significance is 
defined as a/— 21n(L /L max ), where L and L max are the likelihood values returned by the 
fit with the signal yield fixed to zero and at its best fit value, respectively. 

The branching fractions are calculated by using the formula 

^signal 

where N S i gna i, e and N BB - are the number of signal events, the efficiency estimated from MC 
simulation and the number of BB pairs, respectively. We obtain B(B + — > pAtt + 7t~) = 
(5.92lg;ff ± 0.69) x 1(T 6 , B(B + -> pAp°) = (4.78+£g ± 0.60) x 10~ 6 , and B(B + 
pA/ 2 (1270)) = (2.03+°;^ ± 0.27) x 10" 6 with efficiencies of 4.32%, 4.17% and 2.94%, where 
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(a) (b) 




FIG. 2: Distributions of (a) AE (with 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 ), (b) M bc (with 
|A£| < 0.05 GeV) and (c) M n+n ~ (with 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and \AE\ < 0.05 GeV). 
The solid curve represents the fit projection, which is the sum of signal and background (dashed 
curve) estimates. The shaded area represents the sum of signal components i.e. B + — > pA7r + 7r~, 
B + -> pKp° and B + -> pA/ 2 (1270). 



sub-decay branching fractions are included. To verify the validity of the branching frac- 
tions, we select B + -> pA~n + events with 2.10 GeV/c 2 < M A?r + < 2.32 GeV/c 2 for A+ 
and use a 2D (M bc — A£?) extended likelihood fit method. We find the branching fraction 
of B + — > pA c 7r + to be (2.44 ± 0.30(stat.)) x 10~ 6 , which agrees with the value calculated 
using world averages [5j, B(B + -> j9A c "vr+, A+ -> Att+) = (2.25 ± 0.87) x 10~ 6 . 
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TABLE II: Summary. 



Mode 


Yield 


Efficiency (%) 


£(1(T 6 ) 


Significance 


pAir + -K~ 
pAp° 
pA/ 2 (1270) 


167.8l|;° 
131.2±lf;| 

oq -,+14.9 


4.32 
4.17 
2.94 


5.92lg;|| ± 0.69 
4.78±g;64 ± 0.60 
2.03±8;t2 ± 0- 27 


9.1 
9.5 
3.0 



V. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are estimated using high-statistics control data samples. The 
tracking efficiencies are measured with fully and partially reconstructed D* samples. For 
proton identification, we use aA^ prr~ sample, while for K/ir identification we use a D* + — > 
D°7i + , D° — > K~ir + sample. The average efficiency difference for particle identification 
(PID) between data and MC has been corrected to obtain the final branching fraction 
measurements. The corrections are 14.8%, 13.8% and 14.2% for B + — > pA-K + -K~ , B + — > pAp° 
and B + — > pA/2(1270), respectively. The uncertainties associated with the PID corrections 
are estimated to be 4.2% for two protons (one from A decay) and 1-2% for two charged 
pions. 

For A reconstruction, we have additional uncertainties of 3.0%, 2.4% and 2.6% for 
pAir + ir~, pAp° and pA/2(1270), for the efficiencies of tracks displaced from the interac- 
tion point. These uncertainties are determined from the difference between A proper time 
distributions for data and MC simulation. There are also uncertainties of 1.2%, 1.3% and 
1.3% due to the A mass selection and 0.7%, 0.5% and 0.6% due to the A vertex selection, 
for the pAir + ir~, pAp° and pA/ 2 (1270) modes, respectively. 

The uncertainty in the likelihood ratio requirement, TZ, is estimated from a control sample 

it it 7r + ), which has the same number of final state particles. For 
the multiple counting uncertainty, we take the difference of branching fractions with and 
without the best candidate selection as the systematic error. 

Fitting uncertainty is checked by considering different PDF models. MC statistics con- 
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tribute a 1% systematic error. According to our MC simulation study, the charmless B decay 
that most affects our signal determination is B + — > pX7r + 7r~. This mode contributes 1.5%, 
0.5% and 0.5% uncertainties for B + —>■ pAtt + tt~ , B + — > pAp° and B + —> pA/2(1270), respec- 
tively. In the M 7T + 7T - spectrum, we use a threshold function and a phase space shape to model 
the generic B + — > pAn^n^ process. The modeling of the threshold function shape in the 
M^+t,.- spectrum contributes systematic errors of 6.9%, 8.8% and 9.0% for B + — > pAn + n^ , 
B + — > pAp° and B + — > pA/2(1270). Another uncertainty is due to MC modeling of the 
M p -£ distribution and is included by varying the fit parameter of the threshold function and 
adding a phase space component in the fit. The uncertainties due to M p -^ modeling are 1.3%, 
1.5% and 3.2% for B + -> pAn + n-, B + -> pAp° and B + -> pA/ 2 (1270), respectively. We 
neglect the uncertainty due to interference effects between resonances and 7r + 7r~, assuming 
the non-resonant 7r + 7r~ system is predominately in an S-wave. 

We assume that the branching fractions of T(4S) to neutral and charged BB pairs are 
equal. The correlated errors are added linearly and the uncorrelated ones are added in 
quadrature. The total systematic uncertainties are 11.7%, 12.6% and 13.4% for the pA7r + 7r~, 
pAp° and pA/ 2 (1270) modes, respectively, and the details of each decay channel are sum- 
marized in Table [B 



VI. SUMMARY 

Using 657 xlO 6 BB events, we observe non-resonant B + — > pAn + 7r~ and B + — > pAp° 
decays with significances of 9.1 and 9.5 standard deviations as shown in Table HT1 A hint of 
a B + pA/ 2 (1270) signal is found. This is the first observation of a four-body charmless 
baryonic B decay. We also observe low- mass pA threshold enhancements in both the pA^n^ 
and pAp° modes. The branching fraction of the four-body decay is comparable to the 
corresponding three-body decays modes; however, the central values may indicate that the 
hierarchy established in the charmed case still holds for charmless baryonic decay. The 
observed branching fraction for B + — >• pAp° is comparable to that for B° — > pAw~, and is 
about an order of magnitude larger than the theoretical prediction 8J]. 
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